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SUMMARY 

Our Sensitivity-Based Minimal Residual (SBMR) method which is based on our earlier Distributed 
Minimal Residual (DMR) method allows each component of the solution vector in a system of 
equations to have its own convergence speed. Our global SBMR method was found to consistently 
outperform the DMR method while requiring considerably less computer memory. Recently, we 
have developed and tested a new Line SBMR or LSBMR method and a Time-Step-Scaling (TSS) 
method that are even more robust and computationally efficient than our global SBMR method 
especially on highly clustered computational grids in laminar and turbulent flow computations. 

Distributed Minimal Residual (DMR) Method 

This method predicts an optimum amount of correction to each component of the solution vector by 
combining the information from several previous iteration levels. Each of the corrections obtained 
from the past iterations is multiplied by a different weighting factor and these weighting factors are 
determined so that they minimize the overall future residual. Although it is based on general Krylov 
subspace methods, the DMR method [1] differs from them by the fact that weighting factors are 
different from one variable to another in the system. 

Global Sensitivity-Based Minimal Residual (SBMR) Method 

The residual at a grid point depends on the solution vector Q at the neighboring grid points 
including the point itself. The sensitivity of residual (m= 1,2,3: number of equations) with 

respect to (k=typical neighboring points) is That sensitivity can be determined from the 

• 

finite difference equation used in the scheme. For a two-dimensional incompressible flow solved 
using Chorin's artificial compressibility method the solution vector is Q = {p u v) . Suppose we 
have calculated the solution vector Q at iteration levels up to t+n where n is the number of regular 
iteration steps between two iteration levels. Then the change in the solutions between the iteration 
levels can be written as follows. 

Ap k = (pJ + “-(pJ Au^u^-K)*, Av k = KT-(vJ, (1) 
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Using the first two terms in a Taylor series expansion in artificial time direction, each residual for a 
two-dimensional system after n iterations will be 
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Similarly, future residual at t = (t+n)+l can be approximated by 
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Here, a's are the factors that multiply A's to estimate the future solution vector so that it minimizes 
the future residual. Then, (p = (pj + a p Ap k with similar expressions for u and v [23,4,5]. 


Line SBMR (LSBMR) Method : - 

The SBMR method calculates the same a's for the entire computational domain, which cannot 
represent optimum a's for both coarse grid regions and fine grid regions. Therefore, a natural 
conclusion is to allow a's to have different values in the clustered regions. A modified formulation 
(LSBMR) was developed [3,5] to allow these a's to have different values from one grid line to 
another. The LSBMR formulation will be explained using the two-dimensional, incompressible 

flow Navier-Stokes equations as an example. The system has three equations (rmax=3) and three 
unknowns (M=3), that are p, u and v. The acceleration coefficients for those unknowns are ap, a u 

and cty, respectively. If the grid lines are clustered in the j-direction, then each j = constant grid line 

will be assigned its own set of three constant a's. The residual at a grid point (i,j) incorporates a's 

at the neighboring grid points plus the point itself. For the given Navier-Stokes system, this yields 
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On each j = constant line, three values (M=3) of constant a's are determined so as to minimize the 
L-2 norm of the future global residual: 
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for j = 1, 2,- • •, jmax. If jmax is the total number of j = constant grid lines, then substituting the 
equation (4) into (5) results in (jmax) x M algebraic equations for the same number of unknown 
optimum a's. For a given j = constant grid line when using central differencing in j-direction, aJ 

appears only in R^ + (i,j-l), R„ '(i, j) and R^ ,+ (i, j+1). Therefore, we can see that the summation 
over the entire domain leaves the terms only with cJ' 2 , , aJ + ^ and cJ +2 summed along the j 

= constant line. Now we have 
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for j = 12,- • jmax and m = p, u, v. In this example the simultaneous system of equations (6) yields 
a block penta-diagonal matrix equation for (jmax) x 3 optimum a's where each block is a 3 x 3 
matrix. In the general case of a two-dimensional system having M partial differential equations, the 
block penta-di agonal system (6) will have blocks of size M x M. Superior performance of 
LSBMR is demonstrated for laminar (Fig. 1) and turbulent (Fig. 2) flows. 


Time Step Scaling (TSS) Technique 

A simple and efficient additional convergence acceleration technique called the Time Step Scaling 
(TSS) was developed [3,5] for the steady-state solution of the incompressible flow on highly- 
clustered grids. This technique can be considered as an extension of the local time-stepping so that 
time steps can be different for each equation and at each grid point. In the highly-clustered grid 
regions, the grid spacing in the direction of grid clustering is very small and it limits the magnitude 
of Ax to be extremely small for stability reasons (time step is proportional to the square of the grid 
spacing). On the other hand, the Jacobian, J, becomes large in the clustered regions. However, the 
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product JAt is orders of magnitude smaller in the fine grid region than in the coarse grid region. 
Since the transient solutions are not of interest, one may desire to use larger time steps for some 
equations in the highly-clustered region for faster convergence. A close investigation of the 
convergence characteristics of each equation reveals that it is the slow evolution of pressure field in 
the fine-grid region which slows the overall convergence. Therefore, the time step for the continuity 
equation is scaled such that JAt in the clustered regions has the same order of magnitude as that in 
the coarse grid regions. To do this, the maximum value of JAt along the grid line in the direction of 
grid clustering is located and the time step at each grid point along this grid line is scaled as follows 


AtOJ) = c jjj4-j]-( J AtL (7) 

where AT(i, j) is the calculated time step without the TSS, J(i, j) is the Jacobian, c is a user specified 

coefficient (typically 0.8) and (jAx)^ is the maximum value of JAt along the grid line considered. 

For highly-clustered grids, magnitude of JAt in the coarse grid region is orders of magnitude larger 

than that in the fine-grid region. The TSS technique was found to significantly accelerate the 

convergence when highly clustered grids are used (Figures 1 and 2). 

REFERENCES 

1. Lee, S. and Dulikravich, G. S., "Distributed Minimal Residual (DMR) Method for 
Acceleration of Iterative Algorithms," Computer Methods in Applied Mechanics and 
Engineering, Vol. 86, pp. 245-262, 1991. 

2. Choi, K-Y. and Dulikravich, G. S., "Reliability Enhancement of Navier-Stokes Codes Through 
Convergence Enhancement," NASA Propulsion Engineering Research Center Fifth Annual 
Symposium, University Park, PA, September 8-9, 1993, pp. 181-186. 

3. Choi, K.-Y. and Dulikravich, G. S., "Convergence Acceleration of Iterative Algorithms Using a 
Sensitivity-Based Minimum Residual (SBMR) Method," 12th Workshop for CFD 
Applications in Rocket Propulsion, NASA MSFC, Huntsville, AL, April 19-21, 1994. 

4. Choi, K-Y. and Dulikravich, G. S., "Global Sensitivity-Based Minimal Residual (GSBMR) 
Method for Acceleration of Iterative Convergence Rates," Proceedings of the 14th International 
Conference on Numerical Fluid Dynamics, Bangalore, India, July 11-15, 1994. 

5. Choi, K.-Y., "Sensitivity-Based Minimum Residual Methods for Convergence Acceleration of 
Iterative Algorithms", Ph. D. dissertation. Department of Aerospace Engineering, The 
Pennsylvania State University, August 1994. 


180 




Figure 1 Convergence histories for a 2-D straight channel flow including the TSS technique 
(L/H = 5, Re = 1600, 60x60 grid cells, AR max = 10 4 , p = 5, CFL = 2.8, von Neumann = 0.4) 



Figure 2. Convergence histories for a 2-D straight channel flow including the TSS technique 
(L/H = 10, Re = 1.6xl0 6 , 60x120 grid cells, AR max = 10 4 , = 5, CFL = 2.0, von Neumann = 0.4) 
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